Long-range correlations of density in a Bose-Einstein condensate expanding in a 

random potential 
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We study correlations of atomic density in a weakly interacting Bose-Einstein condensate, ex- 
panding diffusively in a random potential. We show that these correlations are long-range and that 
they are strongly enhanced at long times. Density at distant points exhibits negative correlations. 
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The behavior of Bose-Einstein Condensates (BECs) in 
disordered potentials has attracted growing interest of 
physicists during the last few years. In particular, ex- 
pansion of BECs in randonipotentials has been investi- 
gated in detail pj, 13, Isj, Ij, la, l6(. The main interest of 
using BEC systems to study disorder-related phenom- 
ena is that their physical parameters (such as the num- 
ber of atoms, the parameters of the random potential, 
or the strength of inter-atom interactions) can be con- 
trolled quite precisely. During the last two years, consid- 
erable efforts were undertaken to study the atomic den- 
sity n{rA) averaged over disorder in ID [2], 2D [3|], and 
3D [sl, '4| geometries, with special interest in the phe- 
nomenon of Anderson localization [2, 4 Isl, IQ. At the 
same time, little is known about statistical fluctuations 
of ri(r,t) due the randomness of the potential. Recently, 
Henseler and Shapiro ,7| have shown that a BEC ex- 
panding in a random potential is characterized by a com- 
plicated, highly irregular density pattern reminiscent to 
what we know as "speckle" in optics [8[. According to 
Ref. [7| , multiple scattering from the potential completely 
randomizes n(r,t) and reduces the correlation length of 
atomic speckle pattern to a value of the order of the 
healing length of the initial trapped condensate, which 
is the minimal possible correlation length for a coher- 
ent matter wave. The purpose of this Letter is to show 
that the macroscopic coherence of the condensate gives 
rise to genuine interference effects that were ignored in 
Ref. [7|. These interference effects are similar in nature 
to those leading to Anderson localization of the conden- 
sate at large disorder [4l|. They result in stronger fluc- 
tuations and long-range correlations of n(r,i), akin to 
the long-range correlations of conductance fluctuations 
in disordered metals l9\ and the long-range correlations 
of intensity in optical speckle patterns [IC 



11|. For a 



BEC expanding inside an optical waveguide, the long- 
range correlations grow in absolute value with time and 
become dominant in the long-time limit. They can take 
both positive (for relatively close points) and negative 
(for distant points) values. 

We consider a dilute BEC of A^ ^ 1 atoms of mass 
m inside an infinitely long optical waveguide of diam- 



eter d and cross-section A = ttcP/A, parallel to the z 
axis and described by a 2D potential V±{x,y), see Fig. 
[T] The waveguide geometry is rather popular in BEC ex- 
periments and was, in particular, used in the recent work 
on Anderson localization [5J. In the longitudinal direc- 
tion, the condensate was initially confined by a tight ID 
parabolic trap potential Vz{z) = mujlz'^/2 that has been 
switched off to let the condensate expand along the z axis. 
A time T > 1/lOz later the role of interactions between 
atoms in the condensate becomes negligible [2|, |5|, |l3| 
and a weak 3D random potential V{v) is switched on; 
we refer to this moment as i = 0. V{y) is assumed to 
have a white-noise Gaussian statistics: V{y) = and 
V{y)V{y') — u5{y — y'), where the horizontal bar denotes 
averaging over realizations of the random potential. The 
associated mean free path \s t = h'^-K /um? 12] and the 
"weakness" of the random potential is quantified by a 
condition k^H ^ 1, where k^ — ^^/2m^/h and /i is the 
chemical potential of the initial trapped condensate. For 
t > 0, the macroscopic wave function of the condensate 
'0(r, i) obeys the linear Schrodinger equation [l3|: 
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In this Letter we assume that the diameter d of the 
waveguide to which the expansion of the BEC is con- 
strained obeys 27r/fc^ «C d < ^. This corresponds to 
what is known as a "quasi- ID" geometry in the multiple- 
scattering literature [ij]: at distances larger than I, the 

average atomic density f7,(r,t) = |V'(r,i)| is described 
by a ID diffusion process, unlike the recent work [2, |5[ 
where ID (and not 3D) disorder was considered and dif- 
fusive propagation didn't appear. In the quasi-lD geom- 
etry, the localization length at energies of the order of /i 
is typically ^^ ^ ({kf^d)"^, exceeding the mean free path i 
by a factor {k^dY ^ 1. Hence, the condensate expands 
by diffusion until |0| ~ ^^ :^ £ before it starts to be af- 
fected by Anderson localization effects. This regime was 
not accessible in the recent experiment [5| and has not 
been studied theoretically yet. 

For weak disorder k^t ^ 1 and at large distances z':^ (. 
and long times t 3> £/fp (with w^ — hkf^/m being the 




FIG. 1: (color online). Cartoon of a BEC expanding in a 3D 
random potential and confined transversally to a waveguide 
of typical diameter d. We assume 27r/fcp <Si d< I. 
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typical velocity of an atom with kinetic energy /i), the 
average atomic density is independent of x, y and can be 
written as 0, 0, 0] 



n(z, t) 



dk 

2^ 



|0(fc)P 



^e.(2,i), 



(2) 



where |0(fc)|^ (X (l-fcV2fc^)i?(l- |fc| /\/2A;^) @,[l3,[ll, 
with H{x) the Heaviside step function, describes the 
momentum distribution of atoms in the condensate at 
t — Q, tk = h?k^/2m, and P^{z,t) is the proba- 
bility to find a particle of energy e, initially located 
at the origin, in the vicinity of r = (x, y, z) after 
a time t fl2'|. The Fourier transform of the latter 

is Pe{z,n) 



= G 



e+hn/2 



(r,0)G: 



-hn/2 



(r, 0)/27ri/(:, where 



G^{r,r') is the Fourier transform of the Green's func- 
tion of Eq. ^ and i/^ is the local density of states at 
the energy e. In the hydrodynamic limit hfl <^ e, P^ 
is a solution of a ID diffusion equation: P^{z,fl) — 
exp(- |z| v'-«f^/£'e)/2AV-i^£'e L12]> where 1?^ = W£^/3 
is the diffusion coefficient for atoms at energy e. This 
yields 



n(z, t) 



N 



(3) 



where f{x) can be expressed through special functions 
and f{x) ~ 0.6 for a; 3> 1. 

As the condensate expands, its typical size grows with 
time according to (z^) ~ D^t; profiles of atomic den- 
sity are plotted as functions of z in the inset of Fig. [51 
In the main plot of Fig. [2] we show fi, as a function of 
time. The density reaches a maximum at the "arrival 
time" ^arrival — 2z^/I?^ and decays as 1/Vi at long times 
t > ^arrival. This long-time limit is the most interesting 
regime to which we restrict our analysis from here on. It 
is worthwhile to note that Eq. ([3]) breaks down at very 
long times, when Anderson localization comes into play. 
Indeed, Eq. ([3]) predicts fi ^ for i ^ oo, whereas lo- 
calization will "freeze" n at values of the order of N/A^^ 
(for z < S.^), starting from some localization time tioc- 
The latter can be estimated by comparing N/A£_^ with 



FIG. 2: (color online). Ensemble-averaged atomic density n 
of a BEC expanding in a 3D random potential inside a quasi- 
ID optical waveguide. The main plot shows n as a function 
of time for three different distances z. The dashed line is a 
l/\/t asymptote. The inset shows n as a function of z. 



the long-time limit of Eq. ^ or, equivalently, by equat- 
ing (z^) and ^^. One obtains iioc '^ (£/u^)(fc^d)^, which 
shows that a considerable interval of validity exists for 
Eq. ^ between the typical mean-free time (./v^ and 
^loc > ^/v^n when kfj,d S> 1. 

Let us now study correlations of atomic density in the 
expanding condensate, which is the primary subject of 
this Letter. We define the correlation function of density 
fluctuations (5n(r, t) = n{r, t) — n(r, t) as 



C(r,t;r',0 



Sn{r,t)Sn{r',t') 



n(r, i) X n{r',t') 



(4) 



Using the relation V'(r, t) — J d^riG{r, ri, t)(p{ri), we can 
write the numerator of Eq. ^ as 



Sn{r,t)dn{r',t') 



{2TThy 



'tj dtj e-i('i-<^2)t-i(f3-f4)t' 



K (r, t; r', t'- {r,} , {e,}) 0(ri)0*(r2)0(r3)r (r4),(5) 

where the 6-point kernel K is given by the connected 
part of a product of 4 Green's functions, averaged over 
disorder: 

K - G,^ (r, ri)G*^ (r, v^)G,^ (r', v^)Gt^ (r', r4) 



- G,^(r,ri)G*,(r,r2) x G.^lr', r3)G*Jr', r4). (6) 

The largest contribution to K is obtained by de- 
coupling GG*GG* in the first line of Eq. ([6]) as if G 
were a circular complex Gaussian random field: Ki — 
G,^(r,ri)G*^(r',r4) x G,3(r', r3)G*Jr,r2). When K = 
Ki is inserted into Eqs. (O and dUjthe short-range corre- 
lation function Gi studied in Ref. [7| is obtained. Gi is of 
order 1 for r = r', but rapidly decays to zero already for 




FIG. 3: Diagram for the kernel K2 generating the long-range 
correlation function of density fluctuations C2 through Eq. 
([5]). The solid and dashed lines represent the Green's func- 
tions G and it complex conjugates G* , respectively. The 
two parallel G-lines connected by dotted "ladders" symbol- 
ize averages of products of two Green's functionSj_ GG* . 



H 



{tmyZTiK 



ei 



^M- 



Integration should be performed over rs = re = rr = rg. 



|r — r'l ^ 1/fcp. The long-range part of the correlation 
function — C2 — can be obtained by using a next-order 
contribution to K: K2 given by the diagram depicted in 
Fig. [21 This diagram represents an interference process 
between four matter waves that propagate in pairs to 
some point, where they interchange partners before con- 
tinuing to the measurement points r and r'. A proper 
treatment of such a crossing of wave paths is guaranteed 
by the use of the so-called Hikami box diagram [15, 16^ — 
the H box in Fig. [3] — that ensures conservation of parti- 
cle number in the final result. A simplified version of the 
diagram of Fig.[3l corresponding to equal energies tj = ep 
and identical initial points Vj = rg, yields the long-range 
correlation function of intensity fluctuations for a scalar 
wave emitted by a point source in a disordered medium 
llj. Evaluating the diagram of Fig. [3] and inserting the 
result into Eq. ([5]) gives 



— - — -— 2iilh^A 

on(r,t)dn(r,t') — - — „,_ , . x 



dzi 



3m2(27r)4 
2 



j=i 



(7) 

where e± = [e^^ -I- Ck^ ± h{ni — VL2)/2]/2 and ri± = 
±(e/ci — Efcs) + (^1 + ^2)72. This equation can be eval- 
uated numerically and allows for analytical analysis in 
some special cases, as we now show. 

The simplest quantity that Eq. ([7|) allows us to study is 
the correction to the (normalized) variance of the atomic 
density fluctuations: 



(5n2(r, t) 



1 + C2(r,t;r,i), 



(8) 
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FIG. 4: (color online). Position dependence of the correc- 
tion c to the variance of atomic density in an expanding BEG 
for three different times r = tfj,/h. The dimensionless dis- 
tance is ^ = z^k^£/£. The dashed lines show analytic re- 
sults Br^^-^[1 - CiC/r^^'^f] (small () and Dt/C (large C) for 
T = 3 X 10'*. The numerical constants B, C and D given in 
the text. The inset shows the time dependence of c for a given 
(small) C = 0.4. 



where the unity on the right-hand side originates from 
the Ci term and 



C2(r,i;r,<) = 



-AC^^) 



(9) 



with a combinatorial factor 2 added. Here we introduced 
dimensionless variables C — ^\fk^l i and r = i/x/ft that 
feature natural spatial and temporal scales of the prob- 
lem. The dependence of Ci on position <^ and time 
T appears to be given by a universal function c(C, r) 
that does not depend on any parameters of the prob- 
lem. We plot this function in Fig. [3] for three fixed 
values of r. It is quadratic in C, for small C ^ r^/^: 
c(C,r) ~ Bt'^I\\ - C(C/t1/4)2] and decays as DrjC, for 
C > ri/4. Here B ~ 6.3, C ~ 0.3 and D ~ 1.9 are 
constants that had to be calculated numerically. c(C, r) 
grows with time as r"^/^ for ( ^ r^/** as we also show 
in the inset of Fig. |31 This amplification of c{(,t) with 
time can make C2 significant for large t, despite the small 
prcfactor in front of c in Eq. ^. C2(r,t;r, i) becomes 
of order 1 for i - (£/f^)(fc^d)8/3(fc^£)-i/3^ ^hich is stiU 
smaller than the localization time iioc that limits the va- 
lidity of our analysis. 

The long-range correlation of the fluctuations of atomic 
density is obtained from Eq. ([7|) with Az — \z ~ z'\ » 
£. Note that Ci correlation is negligible for such large 
spatial separations. In this work we restrict ourselves 
to equal-time correlations; t = t' . As an example, we 
consider correlations of density at two points symmetric 
with respect to the origin: z = —z' = Az/2, but results 
are qualitatively similar for any sufficiently distant z and 
z' of opposite sign. We obtain 

C2(r, t- r', t) = _-L=^c2(AC, r). (10) 
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FIG. 5: (color online). Long-range correlation function C2 
versus dimensionless distance A(" = 2(" (solid lines) for three 
different times r = tji/h. The dashed lines show asymptotes 
Sr^/'*(1 - BAC/r^'"') and -Ft/A( for t = 3 x lO". The 
numerical constants B, E and F are given in the text. The 
inset shows the time dependence of the correlation for A(^ = 
0.8. 



and nontrivial spatial structure than just a random ar- 
rangement of small regions of high and low density put 
forward in Ref . [7| . Although in the present paper we con- 
sider an uncorrelated, white-noise potential, our results 
can be generalized to correlated potentials in a standard 
way [13, ll8|. We estimate that the long-range corre- 
lations that we study in this Letter should be directly 
observable under conditions of the experiment of Ref. 
[5| for an isotropic random potential with the correla- 
tion length equal to the longitudinal correlation length 
of Ref. [5| and a slightly weaker transverse confinement. 
It would be extremely interesting to extend our calcula- 
tion to the Anderson localization regime (z > ^p), where, 
by analogy with microwave experiments [l9| , one expects 
anomalously large density fluctuations. 

We thank B.A. van Tiggelen and A. Minguzzi for dis- 
cussions. S.E.S. acknowledges financial support of the 
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We show C2 in Fig. [5] as a function of AC for three dif- 
ferent times. For small distances A^ ^ t^^^, the corre- 
lation decays linearly with A^: C2{A(,t) ~ Bt^^^{1 — 
EAC/t'^^'^), where ^ ~ 2. At AC -- r^/-*, it becomes 
negative and reaches a minimum. For longer distances 
AC 3> T^'^, C2 remains negative whereas its magnitude 
decays only algebraically: C2 — —Ft/ AC, with F ~ 1.2. 
In addition to having long range is space, C2 correlation 
grows in magnitude with time (see the inset of Fig. [S]), 
similarly to the variance of 5n. This should facilitate its 
experimental observation. 

Negative correlations of atomic density in an expand- 
ing BEC could be anticipated from the conservation of 
atom number TV which implies J (PrSn{r,t)Sn{r' ,t) = 
and thus requires that the integrand must change sign. 
(Note that our Eq. ([7]) obeys this condition exactly.) The 
important result of our work is to show that negative cor- 
relations occur at large distances between points r and 
r' and that they become increasingly important as the 
condensate expands. Negative correlations of similar ori- 
gin were predicted to exist in reflection of waves from 
a thick disordered slab 17[ . The analysis of long-range 



correlations introduces a new characteristic length scale 
C* - t1/4 or z* - {Df,t ■ i/k^Y/"^ which is much smaller 
than the root mean square size of the condensate (z^)^/^. 
It determines the typical separation between two points, 
situated symmetrically with respect to the initial location 
of the condensate, at which density correlations change 
sign. 

In conclusion, the macroscopic coherence of the con- 
densate prevents the breakdown of the atomic speckle 
pattern ri(r,i) into small independent spots, imposing 
correlations between distant points. As a consequence, 
atomic speckles appear to have a much more complex 
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